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We study a U(N) x U(N) symmetric scalar field model in four and three dimensions. First, using our data 
in four dimensions in the weak coupling region, we demonstrate explicitly that the observed first order phase 
transition is induced by quantum fluctuations. Next, based on the renormalization group and our new simulation 
results in three dimensions we argue that even if the Ua(1) symmetry is restored below the critical temperature 
the QCD finite temperature chiral phase transition for two flavor could be extremely weak first order. 



1. The Model 

Let us consider a scalar field model with action 
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where </> is a N x N complex matrix and D is 
space-time dimension. It has U (N) x U(N) global 
symmetry under transformation: cf> — > U<f>V\ 
where U and V are N x N unitary matrices. 

In contrast to a O(N) scalar model, the sym- 
metry breaking phase transition (PT) for Eq. (|l]) 
can be first order due to the well known Coleman- 
Weinberg phenomenon A qualitative argu- 
ment for the Coleman- Weinberg phenomenon has 
been developed based on the effective potential 
and the renormalization group (RG) |Q. Essen- 
tially first order PT occurs due to the existence of 
run-away RG trajectories. However, physical ap- 
plications of Eq. (fil) rely on the RG flow obtained 
in perturbation [ 2|i or D = 4 — e expansion extrap- 
olated to e = lp|. It is important to extend these 
arguments in a nonperturbative study. 

2. D=4 

Our first numerical results in 4D are reported in 
Ref. Q , which we refer to for details of the lattice 
simulations and discussions of their relevance to 
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the continuum physics. As an illustration for the 
quantum fluctuation induced first order PT, we 
plot in Fig. 1 a comparison of the one-loop per- 
turbative calculation with MC data in the weak 
coupling region. Here the one-loop effect changes 
the second order PT at the tree level to first or- 
der. First order PT were also found in the strong 
coupling region pj, indicating the absence of a 
stable infrared fixed point there. 

3. D=3 

In three dimensions, Eq. (|l]) can be used 
to describe the QCD finite temperature chiral 
phase transition (ChPT) |J. Here the order of 
PT depends crucially on the number of fermion 
flavor Nf and the pattern of symmetry For 
Nf > 3 ChPT is in general expected to be 
first order §. For Nf = 2, if the axial U A (l) 
anomaly remains effective up to the ChPT tem- 
perature T c h, the symmetry of the system will be 
SU(2) x SU(2) - 0(4) and the transition can be 
second order. However, Ua{^) symmetry could 
be restored below T c ^, the symmetry becomes 
U A {1) x SU(2) x SU(2) and ChPT would be first 
order ||. 

So far there is no numerical evidence for a first 
order ChPT for Nf = 2 @. Does this imply that 
the Ua(1) symmetry is not restored below T c h? 
The answer can not be conclusive. In fact, as 
we will show in the following using the effective 
action in Eq. (0): even if we assume the Ua{^-) 
symmetry is restored below T c h, ChPT can be 
very weakly first order. 
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Figure 1. Comparison of the numerical results 
with the bare perturbation calculation (dotted 
line) in 4D. The numerical data are obtained on 
the 4 4 (open circle), 6 4 (open square) and 10 4 
(open triangle) lattices. The solid lines connect 
6 4 and 10 4 data points to indicate the hysteresis 
effects. 



For N = 2 the RG flow for Eq. (@) is plotted 
in Fig. 2 where the /3-functions are obtained in 
D = 4 — e expansion with e = 1. There is a UV 
fixed point at Ai = A2 = (solid circle) and a 
unstable IR fixed point at Ai = 7r 2 e/8, A2 = 
(solid square). The dotted lines in Fig. 2 indi- 
cate the "stability line" (SL). Qualitatively, for 
a renormalized theory located to the left side of 
SL, PT will be first order 0. However, PT can 
be weakly or strongly first order depending on the 
choice of the bare Lagrangian. The RG running 
of Ai , A2 is very slow in the weak coupling region 
and becomes fast in the strong A2 region. A weak 
coupling bare action chosen to the right side of SL 
will need many decades of running in energy scale 
in order to cross SL. Therefore, even though the 
PT will be first order, the jump in the vacuum ex- 
pectation value (VEV), v, at the transition point 
can be many orders of magnitude smaller than 
the "cut-off" scale A and the transition can be 
very weakly first order [Q . 

As an example, we plot 2v 2 — (tr(f>'4>) in Fig. 3 
for Ai = A2 = 0.0625. Obviously, it is almost im- 
possible to decide numerically at this point if the 
PT is second order or weakly first order. Here we 



Figure 2. RG flow in 3D. The arrows point to the 
infrared direction. 



wish to correct an error appeared in the literature 
||. The point at m 2 = —6.8 (solid circle in Fig. 
3) which was claimed to be the critical point of a 
first order PT || is actually deeply in the broken 
phase. 

The search for a first order PT is further com- 
plicated on a finite lattice. The running of Ai, A2 
is controlled by the correlation length £, which 
is limited by the box size £ < L. Therefore, a 
bare action fixed to the right of SL will only be 
able to run across SL in a large enough box 
The first order nature of the PT will not appear 
until the simulation is performed on a large lat- 
tice. There is an interesting example of this phe- 
nomenon in the literature Q|. In the simulation 
of U(N) x U(N) nonlinear sigma model, which is 
in the same universality class as Eq. (|l|), it was 
found that for N = 2 the first order nature of the 
PT could be exposed only on lattices larger than 
f4 3 . We should also point out that for the same 
initial values of Ai , A2 , the RG running is faster 
with increasing N. This explains the stronger 
first order PT for N = 3 in Ref. @. 

To observe the first order PT in a small box, 
one either has to move the bare action close to 
SL in the weak coupling region, or move to the 
strong A2 region where the speed of running is 
fast We give an example in Fig. 4 where the 
bare action is set to be very close to SL. We are 
able to observe clear signal for a first order PT 
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Figure 3. In 3D, (tr(tf(f>)) as a function of m 
at Ai = A2 = 0.0625. The bare perturbation 
prediction is given by the solid line. 



Figure 4. In 3D, Monte Carlo time history for 
tr{^4>) at Ai = -0.22, A 2 = 0.5, m 2 = -1.038. 
Because of the closeness to the "stability line" , we 
are able to observe tunneling between two states 
on a 6 3 lattice. 



on the 6 3 lattice. 

What are the implications for finite tempera- 
ture QCD? Let us assume that Eq. (0) becomes 
a good effective theory for QCD at energy scale 
A < 2itAqcd- Then the strength of first order 
PT crucially depends on the values of the effec- 
tive coupling at scale A. As discussed above, if 
Ai(A),A 2 (A) fall into the weak coupling region, 
finding the first order PT will become a tough job. 
While it is not possible to determine Ai(A), A 2 (A) 
for QCD, it was found in a large color (N c ) ap- 
proximation for the Nambu-Jona-Lasinio (NJL) 
model that Ai = 0,A 2 = 2n 2 /N c \n(A NJL /A) [§ 
where Ajyjtis the cut-off for NJL model and A is 
the scale when the linear sigma model becomes 
a good effective theory. Although NJL model 
is a poor approximation of QCD, we might be- 
lieve that it is at least qualitatively correct. Then 
A 2 can be strong or weak depending on the ra- 
tio A NJL /A (one might take A NJL ~ 4irA QCD ). 
Since A 2 does not depend on Nf , the only feature 
that distinguishes flavor is the speed of RG flow. 
We could immediately draw a conclusion that it 
would be easier to observe the first order PT for 
Nf — 4 than Nf = 2, which agrees with well 
known simulation results g. Even in the case 
of 4itAqcd/ A ~ 1 and A 2 is very large (Eq. (|l]) 
becomes effectively the nonlinear sigma model), 
results of Ref. fij indicate that one may need a 



quite large lattice to expose the first order nature 
of PT. 

Beyond large N c limit, one in principle can not 
eliminate the possibility of a stronger Nf depen- 
dence for QCD. Then for N f = 2, A 2 (A) could fall 
in the very weak coupling region such that PT is 
extremely weak first order. 

I thank R. Pisarski for useful discussions. 
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MC sweeps 



